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Abstract 

We study Chern-Simons (CS) superconductivity in the presence of uni- 
form external magnetic field of arbitrary strength for a system of fermions 
in two spatial dimensions, which are minimally coupled both to the CS and 
Maxwell gauge fields. We have carried out the computation within the mean 
field ansatz. Analysing only the mean field (i.e., ignoring the fluctuations of 
the gauge fields), we find that chemical potential, susceptibility and magneti- 
zation show discontinuities for integer number of filled Landau levels. Taking 
into account the fluctuations of the gauge fields, we find that the masses of 
the excitations increase with the magnetic field, and that the presence of 
nonlinear magnetic susceptibilities show the absence of any critical or pseudo 
critical magnetic field. Finally, an interesting result is that, unlike ordinary 
superconductors, the system is magnetically asymmetric. 
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I. INTRODUCTION 



The relevance of Chern-Simons (CS) gauge theory of planar phenomena like quantum 
Hall effect is by now well established.0 It has also been recognized that the CS interaction 
can lead to a novel kind of superconductivity - characterized by Parity and Time reversal 
(V, T) violation, no Cooper pair formation, two penetration depths in Meissner effect, and 
finally, an antisymmetric (super) conductivity tensor. Proposed originally as a model for 
high-T c superconductors, it attracts continued theoretical interest, partly due to the novelty 
of the mechanism and partly due to the distinct possibility of the existence of such real 
systems. 

This paper is devoted to a study of CS superconductivity (CSS) in uniform finite (ex- 
ternal) magnetic field. Recall that the conventional superconducting phase gets destroyed 
beyond both a critical temperature and a critical field. While there have been extensive 
studies of CSS at finite temperatures (T), there is not much work at finite magnetic field 
(B). We intend to fill this gap here. 

CSS was first established at T = in the pioneering works of Laughlini and Chen, 
Wilczek, Witten and Halperin,! who considered spinless fermions. This was followed by 
an extension to spin 1/2 by Hosotani and ChakravartyS and Chakraborty, Ramaswamy 
and Ravishankar.! The latter found a unique possibility for the existence of CSS with a 
ferromagnetic ground state. There exists also an extensive literature on the T^O properties 
of CSS. See Ref. 6 for details and for references to earlier works. To put it concisely, it was 
found that CSS gives a normal insulating nonmagnetic state beyond a certain temperature. 
The transition to the normal state is over a rather narrow range of temperatures, but does 
not appear sharp enough to qualify unambiguously to be a phase transition. Thus it is not 
clear whether we have at hand a critical or a pseudo critical temperature. 

Based on a mean field (MF) analysis, Hetric, Hosotani and LeJ conclude that there is 
also a pseudo critical magnetic field beyond which CSS would not survive. We propose to 
study this in detail by computing explicitly the one loop effective action for the system in a 



background magnetic field. Further there is the interesting question of the system's response 
when the sign of B is flipped. We anticipate the system to be magnetically asymmetric 
around B = 0. The interesting region around B = —b, (b being the mean CS magnetic 
field), where the particles exhibit net zero mean field will also be examined here. 

Finally, we remark here that the study of CSS here and elsewhere relies on the MF ansatz 
plus perturbative one loop correction. The MF ansatz requires justification, which has been 
attempted in Ref. 6. Alternatively, one might appeal to the success of the MF picture in 
fractional quantum Hall effect (FQHE),! which is in agreement with experiment.! We return 
to a discussion of its validity here as well, although very briefly and only contextually. 

The paper is organized as follows. Sec. IIA displays briefly the formalism of the MF 
ansatz. Mean field results are presented in Sec. IIB. In Sec. IIC, the form factors are 
evaluated. The effective Lagrangian for the magnetic field is obtained by the integration over 
the fluctuating part of the gauge fields in sec. IIIA. The behaviour of low-lying excitations 
are then discussed in Sec. IIIB. In sec. IV, we compute nonlinear magnetic susceptibilities 
and conclude the paper in Sec. V. 

II. MEAN FIELD THEORY 
A. Formalism 

Consider a system of non-relativistic spinless fermions in (2+1) dimensions whose dy- 
namics is governed by the Lagrangian 

C = —F^F^ - ^e^a^ax + ^iD^ - \D k ^\ 2 + - eA p , (2.1) 

where A(a) denotes the Maxwell(CS) gauge field and the covariant derivative = <9 M — 
ie{A ll + a M ). The condition of a fixed density of the fermions is implemented by introducing 
the chemical potential \i (note that we use fi as a space time index as well; this should 
cause no confusion). The last term represents the background neutralizing 'classical' charge 
density. 
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We employ the path integral formalism to evaluate the partition function 

Z = ( [dA] [da] [d^] [di/S] e*/ d " xC . (2.2) 



We proceed with the evaluation of Z by performing the integration over fermionic field first. 
This gives the effective action for the gauge fields incorporating the accumulated effect of 
fermions on the system. The standard method of evaluation of Z is the MF ansatz in which 
one smears out the CS magnetic field to obtain a uniform background (in which the particles 
move). At T = 0, this approach can be justified for large N (N being related to the CS 
coefficient v = Ne 2 /2ir) and for the parameters used in our analysis.!!! In this case the 
external magnetic field will have to also be included along with the mean CS magnetic field. 
The effect of the external magnetic field would depend on its direction relative to the CS 
field. We expand the gauge fields around the configuration 

A = A 2 = ; a = a 2 = ; A x = A^ x = -Bx 2 ; oi = -bx 2 . (2.3) 

Keeping terms up to second order in fluctuations, we find 

Z = Z MF I [dA][da]e iS , (2.4) 



where the mean field action 

%F = ~ % ln Z MF = -* Tr ln (ida-H + ri-lf d 3 x B 2 , (2.5) 
with H = —D 2 /2m. Also, the one-loop effective action is given by 
S= f d 3 x (-\f, v F^ - ^\d„a x 



d 3 x / d 3 y {Ap(x) + a^x)) U^{x, y) {A u {y) + a v {y)) , (2.6) 



where we have represented the fluctuating fields by a and A again. The current correlators 
are given by 



The fermionic currents are given by 

j (x) = e^V , (2.8a) 
3k(x) = -i^ (tfDrf - D^ty) • (2.8b) 

The single particle Greens function G(x,y) = —i(T ip(x)ip^ (y)) can be obtained by 
solving the differential equation 

(id -H + n) G(x, y) = 5®(x - y) , (2.9) 

subject to appropriate boundary conditions. The boundary conditions which we use for 
evaluating G(x , y) will be discussed in the next subsection. T represents the time ordering 
of two fermionic fields. Thus using a suitable limiting procedure one can express fermionic 
current and current correlator respectively in terms of G(x , y) as follows: 

(j (x))=tG(x,x')\ x>=x ^ =t+o+ , (2.10a) 
= tLP" ~ D 'C)G(x,x') xi=x t/=t+Q+ , (2.10b) 



2m 



H-oo(x, y) = ie G(x, y)G(y, x) , (2.11a) 
n ofc (x, y) = ^ [G(x, y)DtG(y, x) - (D y k *G(x, y)) G(y, x)} , (2.11b) 
n„(x, y) = -i^ [DtG{x, y)D?G(y, x) - (D k Df*G{x, y)) G(y, x) 



+DfG(x, y)D* k *G(y, x) - G(x, y)D* k * D?G(y, x) 



-i^-5 kl (5(x -y) + 5{x' - y)) G(x, x')\ x , =x >tl=t+Q+ • (2.11c) 
In the MF approximation the current correlators will be obtained in terms of the Greens 



functions satisfying ( |2.9| ) with the MF configuration. Thus, we pause to discuss the MF 
ground state first before discussing the fluctuations of the gauge fields. 



B. Mean Field Results 



The mean field configuration in this case is rather involved since the external magnetic 
field changes the degeneracy as well as the cyclotron frequency. Since the levels are otherwise 
completely filled, the highest Landau level (LL) is now only partially filled which makes the 
ground state degenerate. 

To handle this, we shall introduce a fictitious spin-like internal degree of freedom, repre- 
sented by the operator U in the Hamiltonian, to the particle which couples to the magnetic 
field and splits the degeneracy. (Alternatively, one could also introduce a background har- 
monic oscillator potential to split the degeneracy. However, we find the former choice more 
convenient). Indeed, if the degeneracy per unit area is p\ = l/2nl 2 (I = [e\b + 5|] -1 ^ 2 being 
the magnetic length of the system), then the 'pseudospin' operator U belongs to that repre- 
sentation which has exactly as many eigenvalues as piA, where A is the area of the system. 
Thus, the Hamiltonian 

H -> H' = H + \lu c U , (2.12) 

where u c = 2npi/m is the cyclotron frequency and A is the dimensionless strength. As 
mentioned, U has eigenvalues given by 

-(^S^^; « i+1 = « < + l. (2.13) 

The modified spectrum for H is shown schematically in Fig. 1. Since A is to be a small 
parameter (which will be switched off at the end of our calculation), it is necessary that 

N B . 

A 1 + x < — ; x = — . (2.14) 
1 1 pA b v ' 

The spectrum of H' is given by 

e ni = + - + Xu^j u c ; n = 0, 1,2, • • • , i = 1, 2, • • • p { A . (2.15) 
The MF Lagrangian (|2.5|) becomes 



1 oo . ji, d2 

£mf = -j E E / t ln ( fco - e ™ + ^ - t ■ (2 - 16) 

n=0 j ~" 

In evaluating the density p = dCy^p/dp, we encounter the integral of the form 

P-. t ~~~^7~7 — 9{ x ) ; (2.17) 

where an additional convergence term exp (ik S) has been inserted. The integral is thus 

{1 for x > 
. Therefore, expression for the density 
for x < 

of fermions follows: 



dC M F 1 



3& 



"l 1 71 n =0 i 



Using the same technique, ko integral in Q2.16 ) can be determined to obtain 



1 00 B 2 
C mf = -j E E - £ni)0{v ~ e„i) - — . (2.19) 

n=0 i 

Note that the chemical potential (which at T = equals the Fermi energy) can be now 
determined. In the limit A — > 0, 



„H J * +1MK , (2.20) 

where K = N/\l + x\. The upper case corresponds to fractional filling factor K and the 

lower case to the integer values of K. [x] denotes the largest integer less than x. Note that fi 

is discontinuous at those values of B, where it corresponds to integral number of fully filled 

LL. In between the two integer filling fractions \i varies linearly with B as u c varies in the 

same way. In fact, u c (x) = u c (0)\l + x\. Fig. 2 shows the variation of the chemical potential 

with the application of external magnetic field parallel to the mean CS magnetic field. 

The MF Lagrangian for K fractionally filled levels becomes 

1 ( [K ~ 1] \ B 2 

£mf = -j E E - **) + E (t* - W) - ■ (2-21) 

71 \ n=0 i i<i J Z 

Here iq denotes the quantum number of the highest occupied level corresponding to the 
eigenvalue Ui = Uq of U. The summation over n and i can be easily carried out. In the 
limit A — > 0, ( |2.2lD gets a simple form 



£mfH -i 6+ ™ 2 + ™-- . ( , 22) 



The last term in Eq.( 2.22|) is due to the kinetic term of the Maxwell field and does not 



contribute in the calculation of electro-magnetic response. The MF magnetic susceptibility 
can be readily evaluated (by the omission of Maxwell kinetic term) as 

Xmf = oo9 = < ■ (2.23) 



as 2 



2-Km 



This is also discontinuous when [K] passes from one value to the other. This is the well 
known de Haas-van Alphen effect for diamagnetism. Recall that susceptibility diverges 
quadratically for B — > in the conventional de Haas-van Alphen effect where there is no 
internal field.0 Similarly, here, it diverges for B — > —b. The variation of Xmf with the 
magnetic field is shown in Fig. 3. We note that a similar behaviour was also obtained by 
Hetric et al.@ in their calculation of MF magnetization, which we obtain as 

B 



M MF = -^(l+2[K)-2 [ ^±([K) + l) 



1 b 



sign(& + 5). (2.24) 



However, they do not consider the fluctuations of the gauge fields which should play an 
important role in the response of the system to the applied magnetic field. We shall pursue 
this issue in the next section. 



C. Current Correlation Functions 

The evaluation of the form factors requires the Green function ( |2.9|) . We first define the 
frequency transformed Greens function G U (X , X') to be 

G(x, x') = [ pe-^-^G^X , X') , (2.25) 

J ATI 

which clearly satisfies the differential equation 

[u-H' + fi]G u {X, X') = 5 {2 \X-X') . (2.26) 

o 



This can be solved, as usual, in terms of the complete set of eigen functions ip n i(X) of H' 
giving 

G(x ,*')=/- E • (2.27) 

ic* 27r ^ w - e ni + /i 

The contour C for the frequency integration has to be chosen so that G(x , x') satisfies the 
boundary conditions 



G(x,x') ~ <{ -° (2.28) 

En>^ £, e-*^-")^) + £ J>io e -<( e i*i*-'0(*-'') , f > f . 

More explicitly, the contour C must pass below the poles at to = e ni — \i for n < [K — 1], for 
all z and n = [K], i < iq and above otherwise. 

We now evaluate current correlation functions in the momentum space following the 
elegant procedure given by Randjbar-Daemi et al.0 Gauge and rotational invariance requires 
that the current correlator has the form 

IP> , (?) = U (u , (?) (q 2 g^ - q»<f) 

+ (II 2 - n ) (u , f)5^5 vj (f5 ij - gV') + ini(w , gV^g A , (2.29) 

where the last term is the parity and time reversal violating contribution. Here q 2 = uj 2 — q 2 . 

For the purpose of evaluating low energy effective Lagrangian, it is sufficient to compute 
the form factors Ho, Hi and H2 in the limit uo — > 0, q 2 — > 0. The limits commute in the 
evaluation of form factors. For simplicity, we take the limit q 2 —>■ first. In the limit q 2 —>■ 0, 
the form factors are given by 

)0) = ^f^ E (n^! + (n + 1)^) ^ ( 3Qa) 

IIi(w, 0) = U u c , (2.30b) 



2-irmpiA J ^ (uj' - e mj + /i)(u/ - e ni + fi-u) 

x n(n - 1)6 

n,m+2 ~\~ 3fl ^n,m+l 

+ (2n + l) 2 ,5 

11 in 

+ 2(n + l) 2 5n,m-i + (n+ l)(n + 2)5 n , m _ 2 j . (2.30c) 



The integral which we encounter in the evaluation of (|2.30| ) is 



du' 



1 



for 



2ni (u)' - e mj + ji) (u 1 - e ni + /j, - to) 

m > [K + 1] , for all j ; m = [K] , j > i ; 
n < [K] , for alH ; n — [K] , i < i 

n > [K + 1] , for all i ; n = [K] , i > i ; 

-(eni ~ tmj + ^) ; for < 

m < [K] , for all j ; m = [if] , j < io 
otherwise . 



(2.31) 



Thus, the form factors (in the limit u — > , — > 0,) are obtained as 

e 2 K 2 



e 2 K ^ e 2 K 

n = - — , n x 



2-KUJr 



— , n 2 

2% 2nm 



(2.32) 



Note that they are dependent on the applied magnetic field. In terms of their values at 
B = 0, they are expressed as follows: 

e 2 mN 2 n (0) 



47T z p |1 + X\ \1 + X\ 



ni(x) 
n 2 (x) 



e 2 N 



Hi(0) 



2/T |1 + x| |1 + x| 
e 2 AT 2 n 2 (0) 



27rm |l+x| 2 |l + x| 2 



(2.33a) 
(2.33b) 
(2.33c) 



The behaviour of the form factors for negative x, specially for x = — 1, will be discussed 
later in section (IV). 



III. GAUGE FIELD FLUCTUATIONS 

A. Effective Lagrangian 

Now we choose the Coulomb gauge, diAi = 0; didi = in order to evaluate the fluctuating 
part of the partition function ( |2.4j ) by an integration over the gauge fields. Basically one 
computes the contribution to the partition function from the effective gauge field modes 
(collective excitations). We obtain 



ln Z eff = 2 Tr ln 



2 2 W 2 2 



(3.1) 



The collective modes have the dispersion relation u = oj±(q 1 ), with 

1 



2d 



C 2 ± v/C| - AC X C 2 , 



and 



Ci = n 2 , 

c 2 = n (n + n 2 )«f + z/ 2 (n 2 + 2n ) + (y- n ly 
c 3 = noiiag 4 + [z/ 2 (n + n 2 + n n 2 ) + {v - u^ 2 } f + v 2 u\ . 



Therefore, one finds the effective Lagrangian 

i f°° du f d 2 q 



ln 



{u 2 - u 2 + ){u 2 - u 2 _) 



" eff 2 7-oo 2vr J (2tt) 2 
The divergent frequency integral can be regularized as in Eq. ( |2.18|) to obtain 



C 



eff 



d 2 q 



{u + + ojJ) . 



(3.2) 



(3.3a) 
(3.3b) 
(3.3c) 



(3.4) 



(3.5) 



This is the effective Lagrangian for the system coming from the fluctuations of the gauge 
fields, where the dependence on the external magnetic field arises through the form factors in 
the dispersion relations. Below we discuss the behaviour of u± as functions of the magnetic 
field. 



B. Low- lying Excitations 

We first consider, for simplicity, the neutral system (i.e., excluding the internal Maxwell 
gauge field). The dispersion relation that we obtain, for the neutral system, from the 
equation of motion of the CS gauge field is given by 

ITo (IW 2 - Tl 2( f) -(u- nO 2 = . (3.6) 
1 1 



Recall that in the absence of an external magnetic field, the tree level CS term exactly 
canceled with the dynamically generated CS term (i.e., v = Hi). Thus, the system possessed 
a super-fluid mode with the low-lying massless phononic excitation, 

J = . (3.7) 

(This in fact is the collective pseudo-Goldstone mode which is absorbed by the Maxwell 
gauge field upon coupling the particles to the EM field thereby producing a gap and leading 
to the Meissner effect). 

When we introduce external magnetic field, this mode acquires a mass (M), 

M 2 = ^ c 2 (0)x 2 (l + a;) 2 , (3.8) 

as the cancellation of the CS term no longer holds and hence the super-fluidity diminishes. 
The dispersions are shown in Fig. 4 for three different values of B. We have choseniS iV = 10, 
e 2 = 10 5 cm _1 and the value of two dimensionless parameters e 2 /m and p/me 2 to be 10 -5 
and 10 _1 respectively. Observe that the mass gap increases with the value of B although 
the dispersions are parallel. Thus, the velocity of the mode decreases with the increase of 
B. 

We now come back for a discussion of excitations of the charged system. The low-lying 
excitations of the system are two massive photonic modes. These two modes are nothing 
but lu± given in equations (|3.2| and ^3|) . Fig. 5 and Fig. 6 show how the frequencies disperse 
with momentum for three chosen values of B for the modes. The masses M 2 increase with 
B. Note however, that the damping lengths = M± 2 have to be distinguished from the 
penetration depth,!'! which in fact decreases with increasing B, signifying the absence of 
the Meissner effect. This is seen from the calculation of the magnetic susceptibility below. 



IV. NONLINEAR RESPONSE 



Here we study the nonlinear response of the system to an external field by computing 
the higher order susceptibilities defined by 

I rt 



where F is the free energy density of the system which can be computed from the Lagrangian 
( p.5|) as F = — C q q and is recognized to be the linear response susceptibility. Note that 
the complete response of the system is obtained by adding the corresponding MF values 
which may be inferred from Eq. ( [2.23 ). 



We expand the sum u + + u>- up to order (f, since it is a low energy- momentum theory, 
and obtain 



uj + + u;_ = Fi + F 2 q 



->2 



(4.2) 



with 



Fx = v 
F 2 



2 1 / IlA 2 2Ui 
1+ — + ^(1-— ) 



1 



2F 1 



n ' n 2 
, . n 2 



n, 



^ (n + n 2 + n n 2 + (i - ^ 



(4.3a) 
(4.3b) 



We introduce a momentum cut-off A to obtain 



F = — 

8tt 



FxA 2 + -F 2 A 4 
2 



(4.4) 



The susceptibilities x > which are apparently a function of A, are calculated to be 

A 2 



X (r WA) 



d r+2 F 1 A 2 <9 r+2 F 2 



S7r(r + 1 

To determine A, we demand that the linear magnetic susceptibility 

d 2 F 



(4.5) 



X 



(0) 



dB 2 



•1 



£=0 



which we know already and independently from linear response theor 
be 

2p \ 



j7T 2 P 

AT 



2 mNj m 2 N 



(4.6) 

1 21 This fixes A 2 to 
(4.7) 



This value of A is appropriate as a cut-off momentum for high magnetic field also since 
A 2 / 2 »s 1-8/1 1 + x\, using our chosen values of parameters. Therefore, unless x is very high, 



A is quite reasonable as a cut-off momentum. The cut-off independent non-linear magnetic 
susceptibilities can now be readily evaluated from ( [4.5| ) by the substitution of the value of 
A 2 ( |4.7p . They can be analytically determined, to a high degree of approximation to be 

X W = (-l) r+l! ^ • (4-8) 

Notice that higher order susceptibilities are nonvanishing and do not even have numerically 
small values, which clearly shows that there is neither a critical nor a pseudo critical field 
which would characterize the phase of the system. However, as B — > oo, the free energy 
has only a linear dependence on B, which means that all the susceptibilities vanish. The 
system returns to its normal state asymptotically. We also notice that since both even and 
odd order susceptibilities survive, the system is magnetically asymmetric around B = 0, as 
a consequence of the V, T violation inherent in the theory. 

It is instructive to study the behaviour of free energy as a function of applied magnetic 
field. It has the (approximate) form 

b 2 b 2 

F(B) « — — - — —--(1-B/b). (4.9) 
v ; 2(1 + B/b) 2 v ' ' K ' 

Clearly, the leading order contribution comes from B 2 since, as we know from the linear 
response analysis,i@ the system does not possess any spontaneous magnetization. The 
curvature of F(B) at B = gives the linear susceptibility. The behaviour of F(B) is shown 
in Fig. 7. It may be seen (see also Eq. fl4.9Q ) that the free energy diverges at B — —b. This 
requires some discussion, as it corresponds to a zero mean field. 

As B — > —b, the inter Landau level spacings squeeze and hence the energy spectrum 
approaches the continuum; more and more number of LL will be filled up even as the 
particle density remains the same. Therefore, the values of the form factors ( |2.33| ) increase 



and they diverge leading to the divergence in F. The correctness of the above observations 
hinges crucially on the validity of the MF ansatz in this region. Some indirect support for 
the validity can be obtained from the related phenomena of FQHE, where the composite 
fermion model! predicts a zero mean field at the filling fraction v = 1/2, thus leading to a 



free fermion like system. The analysis of Halperin, Lee and and a recent experiment 

by Du et al.i do support this prediction. More specifically for the system at hand, we recall 
the argument given in Ref. 6. One simple criterion is that MF would be plausible if the 
interparticle seperation is less than the magnetic length I which is a measure of the scale 
over which the single particle wave function extends, i.e., 

^p-i </ — — >2tt . (4.10) 

Therefore, the MF theory will be invalid for very high \x\. However for moderate \x\, the MF 
ansatz works very well. Note that Q4.10 ) suggests that the MF theory is valid at x = — 1. 



Before we end this section, we remark that the MF value for F approaches a constant 
value, with an ever increasing number of discontinuities as x — > — 1. In contrast, the contri- 
bution from the fluctuations diverges linearly and thus dominates over the MF contribution. 
To illustrate this, if the generalized susceptibility defined as x( x ) = dM.{x)/dB is con- 
sidered, it is easy to see that the fluctuation part diverges as -1/(1 + xf unlike the MF 
contribution which diverges only quadratically in the same limit. Finally, note that even 
around x — 0, the susceptibilities get a dominant contribution from the fluctuations. 



V. CONCLUSION 

To conclude, the mean field properties such as chemical potential, magnetization and 
magnetic susceptibility oscillate as functions of the external magnetic field. Considering the 
fluctuations of the gauge fields we find that the masses of the mode of excitations increase 
with magnetic field. The nonlinear susceptibilities arising from the fluctuations of the gauge 
fields are computed which show the absence of any critical or pseudo critical field. Further, 
there is a unique asymmetry in the system at B = as far as the magnetic properties are 
concerned, and would be an interesting property to look for in case candidates for CSS are 
proposed in real systems. 
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FIG. 1. The energy levels for MF ground state. N Landau Levels are filled at B = 0. At B ^ 0, 
there are K fractionally filled (topmost level is not fully filled) levels. Split levels arise due to the 
switching on of the fictitious interaction. Each split level can accomodate only one particle. The 
number of split levels for each LL are equal to the degeneracy of each LL. The horizontal dotted 
lines represent the Fermi levels. 

FIG. 2. Chemical potential \i (solid line) is plotted against the applied magnetic field B, N = 6, 
in the units of 2'Kp/m = /i(B = 0). B is in the units of CS magnetic field b. Note that fx shows 
discontinuities at those values of B where the number of filled levels are exactly integral in number. 
Also, note that the values of fi at integer filling (denoted by dark dots) are the same and equal 
that at B = 0. 

FIG. 3. MF magnetic susceptibility XMF ( son d lines) is shown, for N = 6, against B/b. It also 
shows the discontinuities in those values of B for which filled levels are integers. Its values for 
those values of B are shown by dark dots. 

FIG. 4. The dispersions of the phononic modes for different values of B are shown. Here 
lo = | et> | 1 / 2 is the magnetic length of the system in absence of B. The numbers associated 
to each curve are the applied magnetic field B in units of b. Note that the mass gap of the mode 
increases with B. 

FIG. 5. The dispersion relation uo\ as a function of q 1 for different values of B. The numbers 
associated to each curve are the applied magnetic field B in units of b. 

FIG. 6. The dispersion relation uo 2 _ for different values of B. The numbers associated to each 
curve are the applied magnetic field B in units of b. 

FIG. 7. Free energy is shown as a function of magnetic field. 
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